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Abstract 

We evaluate analytically higher terms of the e-expansion of the three- 
loop master integrals corresponding to three-loop quark and gluon form 
factors and to the three-loop master integrals contributing to the elec- 
tron g— 2 in QED up to the transcendentality weight typical to four-loop 
calculations, i.e. eight and seven, respectively. The calculation is based 
on a combination of a method recently suggested by one of the au- 
thors (R.L.) with other techniques: sector decomposition implemented 
in FIESTA, the method of Mellin-Barnes representation, and the PSLQ 
algorithm. 
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1 Introduction 



One year ago a method of multiloop calculations [JJ based on the use of dimensional 
recurrence relations [2J and analytic properties of Feynman integrals as functions of 
the parameter of dimensional regularization, d = 4 — 2e, was suggested (the DRA 
method). Then it was successfully applied [3] to the analytic evaluation of the terms 
of order e° of the e-expansion of two most complicated master integrals corresponding 
to the three-loop quark and gluon form factors. These results provided the possibility 
to obtain completely analytic expressions for the form factors. The power of the DRA 
method was further demonstrated in Ref. [1] where the three previously unknown 
terms of the expansion in e, i.e. up to transcendentality weight ten, of the three-loop 
non-planar massless propagator diagram were evaluated. Some additional details of 
the method were presented in Ref. [5]. 

The main point of Ref. [2] was to emphasize that once analytic expressions, in 
terms of well convergent series, are obtained within the DRA method it is easy to 
evaluate extra terms of the expansion in e. In the present paper we further exploit 
this nice feature of the method by evaluating, up to the transcendentality weight 
typical to four-loop calculations, higher terms of the e-expansion of the three-loop 
master integrals corresponding to three-loop quark and gluon form factors and to the 
three-loop master integrals contributing to the electron g — 2 in QED. 

In the first of these two problems, all the master integrals apart from three most 
complicated ones were evaluated in Refs. [6[ [7J. In fact, the word evaluated means 
here the evaluation up to the order of e which appears in the finite part of the form 
factors. Mathematically, this means the evaluation up to transcendentality weight six. 
Then one of the three most complicated master integrals (called Ag t i in Refs. P, [7] 
and in the present paper) and the pole parts of A 94: and A) j2 (shown in Fig. 1 in 
the next section) were evaluated analytically, while the e° parts of Ag^ and A 9i2 
were evaluated numerically — see Refs. jH |9]. These two missing ingredients were 
calculated in Ref. [3]. Motivated by future four- loop calculations, the authors of 
Ref. [TU] presented one more term of the e-expansion for all the master integrals but 
Ag t i, Ag^ and Ag^. For Ag t i, this was done in [1]. In the present paper we evaluate 
all the master integrals up to transcendentality weight eight which is intrinsic to a 
four-loop evaluation of the form factors. Explicitly, the constants present in results 
for highest powers of e are linear combinations of 7r 8 , Cf 71 " 2 , C5C3 an d C-6,-2- 

The second family of master integrals we are evaluating is connected with the 
anomalous magnetic moment of the muon, i.e. g — 2 factor. The three-loop evaluation 
was performed in Refs. [UJ [12]. The evaluation at the four-loop level, within a 
pure numerical approach where each of the Feynman integrals involved is evaluated 
numerically, without a reduction to master integrals, was already performed — see 
Ref. [13J and references therein. Even the numerical evaluation in five loops has been 
started — see a very recent paper [T5] and references therein. 

In Ref. [15] the status of partially analytic evaluation of the four-loop g — 2 factor 
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was presented. In this approach, an integration by parts (IBP) reduction [16] to 
master integrals is used. However, it was pointed out that the analytic evaluation 
of the corresponding master integrals was not possible for the moment. To evaluate 
the master integrals numerically with a high precision it was supposed to use the 
method developed in Refs. (TTJ HE] and based on difference equations. Keeping in 
mind that all the three-loop g — 2 master integrals will inevitably appear in a four- 
loop calculation because they are present in factorizable four-loop diagrams where 
three-loop master integrals enter in products with one-loop integrals having poles in 
e, Laporta calculated numerically [19] several higher order terms of the e-expansion 
of the three-loop g — 2 master integrals, with the use of the method of Refs. p21 [18] • 
We are more optimistic about the possibility to evaluate four-loop g — 2 mas- 
ter integrals analytically. To make a step in this direction we have performed an 
analytical calculation corresponding to the numerical calculation of Ref. [19]. More 
precisely, we have evaluated analytically higher terms of the e-expansion of the three- 
loop g — 2 master integrals up to transcendentality weight seven which is intrinsic to 
a four-loop evaluation of the electron g — 2. So, the highest terms we present involve 
(7, , Li7(l/2), In 7 2, C— 5,1,1? e ^ c - m our calculation we use the DRA method, obtain 
results for the master integrals with a very high precision (~ 350 — 500 digits) and 
then apply the so-called PSLQ algorithm [20] to arrive at analytical values. There 
are fifty four independent transcendental constants in the four loop basis so that this 
high precision is mandatory. (Typically, one needs the accuracy of at least seven 
digits per constant. Observe that this high accuracy is not provided by the method 

of Ref. pans].) 

The three-loop g — 2 master integrals were analytically evaluated in Refs. [TH [121 
[2H [22] • In fact, in [2H 122], another problem was considered: the evaluation of three- 
loop renormalization constants, for which master integrals are almost the same, up to 
one additional master integral. This three-loop evaluation was typically performed up 
to transcendentality weight five, with some exceptions where some master integrals 
were expanded in e up to weight four. 

In the next two sections we present results for the master integrals for the three- 
loop form factors and for the g — 2 factor, respectively. 

2 Master integrals for quark and gluon 
form factors 

The master integrals for the three-loop form factors are shown in Fig. [TJ Two external 
momenta are on the light cone, p\ = p\ = 0. In Ref. [I] the DRA method was applied 
to the evaluation of these master integrals. The corresponding results for arbitrary 
d are obtained in the form of multiple well-convergent series. Here we present our 
analytical results on the higher terms of the e expansion which were obtained by 
calculating these series with a high precision and then applying the PSLQ algorithm 
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Figure 1: Master integrals for the three-loop form factors. Internal lines denote 
massless propagators l/(k 2 + iO). 



The high orders of e-expansion allowed us to guess, for all master integrals, ex- 
cept A 64 and the three most complicated ones (Ag t x, A 9t2 , and A 9)4 ), the factors 
which makes the expansion homogeneous in the transcendentality weight. Taking 
into account the fact that these four master integrals can be replaced by the integrals 
with numerators [9] (or denominator squared for Aq^), see Fig. [2j which are also 
homogeneous in the transcendentality weight, we have a complete basis of master 
integrals with homogeneous transcendentality weight. Let us present the results for 
the expansion of the integrals in this basis. We arrange the results in the order of 
increasing complexity level, the notion introduced in Ref. [I]. The loop integration 
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Figure 2: Master integrals with a homogeneous transcendentality weight to replace 
the corresponding integrals without numerator. 



measure is taken as 

r\ d h 
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m d/2- y > 

Zero complexity level 

Master integrals with the zero complexity level nullify upon shrinking of any internal 
line. As we have already argued in Ref. pi], they are always expressed in terms of the 
r-functions. We present their expansion here only for convenience of the reader. 
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Non-zero complexity level 

Integrals with complexity level equal to n > are expressed in arbitrary dimension 
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d in terms of an n-fold series. We present these integrals in the same order as in 
Fig. [TJ We also present results for Aq A , Aq^ Ag 2 , and Ag 3 which have the uniform 
transcendentality and can be used as master integrals instead of A^, Ag t i, Ag^, and 
Aq3, respectively. 
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Here 7 is the Euler constant, ( m = C( m ) ; an d Cmi,m 2 = C( m i? m 2) are multiple zeta 
values (see, e.g., Ref. [23] ) 

00 2i — 1 ifc-i — 1 fc / 
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3 Master integrals for the g — 2 factor 

The master integrals for the three-loop g — 2 factor are shown in Fig. [3j The solid 
internal lines correspond to massive propagators l/(—k 2 + m 2 ), while the dotted lines 
correspond to massless ones, l/(— k 2 ). The external solid lines denote the incoming 
and outgoing on-shell momentum p, so that p 2 = m 2 . We set m = 1 for simplicity. 
Again, the ordering in Fig. [3] corresponds to increasing complexity level. 

It is known [TT1 [T2] that at the three-loop level the g — 2 factor is determined by 
seventeen master integrals. In fact, in an earlier paper [TT] there was one more master 
integral (In in the notation of [H]) having the same set of denominators as - 
see Fig. HJ In a later paper [12] it was recognized that this additional integral can 
be expressed via G 44 and G 3 , at least, up to a high power of e. The DRA method is 
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Figure 3: Master integrals for the three- loop g — 2 factor. Each line corresponds to a 
definite complexity level indicated to the left. 




hi 



Figure 4: The integral In (left) and the auxiliary diagram used for its reduction 
(right). 

based on the consideration of relations exact in d, so that we needed an exact relation 
between In, G 4i4 , and G 3 . For this purpose, we have performed an IBP reduction of 
the integrals with denominators determined by the diagram depicted in Fig. HI This 
reduction gave us the expected relation: 

Id — Pi 1 
hi(d) = ^-^CAd) - -G 3 (d) . (26) 

Since we are oriented at a future four-loop calculation, we consider two more inte- 
grals, G 5j 4 and Gj, which also appear in the factorizable four-loop master integrals. 
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The application of the DRA method allowed us to express all master integrals in 
d dimensions in terms of one-, two-, and threefold series with maximal nested depth 
equal to the complexity level of the integral. These representations are available upon 
request from the authors. Here we present expansions of the master integrals near 
d — 4 obtained by numerical summation of the series followed by the application of 
the PSLQ algorithm. Again, for convenience of the reader, we start from listing the 
results for zero complexity level master integrals expressible via T-function. 
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+ 8(,3 + 



^ 1 ^> ! + '""( ^ + 32C 3 2 1 + A + !)(irr\;, + 1992C 7 



l-3e)(2-3e)(l-2e)(l-e) 1 e 3 3e ^ 45 

J V 945 ' "'V ' " V 45 
+0 (e 5 ) ) , (29) 



2 1 - 2 T(2-e)r(e-l) 2 r e-i T(3e-2) 
G 44 4 - 2e = M . ; . 2 -!— 

r (2c - |) 

(l-4e)r(l + e)3 ( 2 I6C3 4^ 

" (l-3e)(2-3e)(l-2 e )(l- e )[3e 3 3 15 ^ 5 

iw-T) + K-^ %+i32s(r ) +o ^} (3o) 



— e 
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r(3-4e)r(2-2e) 



2 



(l_ 6e) r(l + 6) 3 J 1 2tt 2 38C 3 151tt 4 6 2 / 767T 2 C 3 , w 
(l-4e)(l-2e) 3 ] 3e 3 3e 3 90 € \ 3 ^ 5 



-e 3 



(1^! + !|C|) _ £l (!^ + 860A5 + 12434C7 ) +0(e »)} (31) 



~ , xt,/ 3 23 105 /275 4tt 2 \ / 567 , 

Gli3 (4-2 £ )=r(l + £ fJ S 3 +i? + - + (- + -j + e(-- + 1 ^ 

- ln2 + 28C,) + e 2 + ^ - 60tt 2 ln2 + 210( 3 - ^ + 16n 2 ln 2 2 

+ 8h«2 + 192 J + e 3^_144015 + 385^ _ ^.^ + ^ _ s£ 
y \ 128 2 3 

124 48 ln 5 2 
+120tt 2 ln 2 2 + 60 ln 4 2 + 1440a 4 + vr 4 ln2 - 32tt 2 ln 3 2 + 1152a 5 

J^fSl _ 930C5 ) + e'f-H5^ l ^ _ „ 55l 2 ln2 + ™^ _ ^ 

3 V V 256 12 2 18 

+580tt 2 ln 2 2 + 290 ln 4 2 + 6960a 4 + 62vr 4 ln2 - 240vr 2 ln 3 2 - 72 ln 5 2 + 8640a 5 

-100vt 2 C 3 - 6975Cs - - — vr 4 ln 2 2 + 48vr 2 ln 4 2 + + 6912a 6 

135 5 5 

, A ,/ 7710087 18585vr 2 8281 2 , 
+2880s 6 + 80vr 2 ln2 C 3 - 1220C| ) + e ( ^2~ + § ^ 

+ 5^7C3 _ 2387^ + 23io?r2 ln22 + ln42 + 27720a4 + 899^ ^ 

-1160. 2 ln 3 2 - 3481n'2 + 41760a 5 - - ™^ - ^ - 186. 4 ln 2 2 

3 2 9 

+360tt 2 ln 4 2 + 72 ln 6 2 + 51840a 6 + 21600s 6 + 600vr 2 ln2 Cs - 9150C 3 2 

784 fi , n 248 4 , o 288 2l , 288 ln 7 2 55680, rt 

+ vr 6 ln2 + 7T 4 ln 3 2 vr 2 ln 5 2 + 41472a 7 ln2 s 6 

45 5 5 35 7 

4300tt 4 C 3 nin oi 2 «> 69600, n , 2 32086vt 2 C 5 2 « > 

+ ^— ^ - 240tt 2 ln 2 2 C 3 + — =— ln2 C 3 + =— ^ + 16740 ln 2 2 ( 5 

Do ( ( 

579651(7 + 55680s 7a _ 65280 V . + Q ^ . (32) 



^) +0(e «)}. 
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3 i 2 23 35 275 / 189 112Ca 



^ 5 (4- 2e ) = r(i + e )^- + - + - + - + e -- + 



-e 2 



( 14917 ™^ ^vr 4 32 2l 2n 321n 4 2 nr „ \ 



o/ 48005 4060C 3 68tt 4 2i 2 , 4 272 , 

-e 3 -— + — — + 80vr 2 ln 2 2 - 80 ln 4 2 - 1920a 4 - — tt 4 ln2 

\ 32 3 3 15 

64 2l „ 641n 5 2 , ro „ 4 / 1108525 r0rtrt> 986tt 4 

__ 7r 2 ln 3 2 + __ ^gg^ + 1240Cs \ _ £ 4/ ___ 53goC3 + ___ 

+il^7r 2 ln 2 2 - 11601n42 - 9280a 4 - 136tt 4 ln2 - 160tt 2 ln 3 2 + 96 ln 5 2 
3 3 

^9-7r 6 272 641n 6 2 
-11520a 5 + 9300Cs + — + — 7T 4 ln 2 2 + 32tt 2 ln 4 2 - - 9216a 6 

5 5 5 

4880C|\ 5 /2570029 57967C 3 , 1309I 4 , 1CJ „ 2 , 2r 



-3840s, + j - S ^ + — + 1540^ hrt 

1972 2320 

-1540 ln 4 2 - 36960a 4 7r 4 ln2 n 2 ln 3 2 +4641n 5 2 - 55680a 5 

3 3 5 

+44950(5 + 48tt 6 + 408tt 4 ln 2 2 + 240vr 2 ln 4 2 - 96 ln 6 2 - 69120a 6 - 28800s 6 
+ 122C0g-^l»2 -^V2 -12^2+^-55296^ 

74240, „ 720tt 4 C 3 92800, rt >2 130360tt 2 C 5 rtrt0rtrtl 2 « > 

+ n2 s 6 — n2 Cl — - 22320 n 2 2 C 5 

7 7 7 3 21 

+ 7 = _ 742^ + 87044 Jf!5 X o(e8) \ ^ 
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1 7 253 2501 / 59437 64tt 2 
G 4 , 6 (4-2 e )=r(l + e) 3 - + — + — + + 



e 3 2e 2 36e 216 V 1296 9 

.,, 2831381 2272tt 2 256 2 , 1792C 3 \ ,/ 117529021 
-e 2 — — + 7T 2 ln2 + — - e 3 ' 



7776 27 3 9 J \ 46656 

49840tt 2 9088 2 , 63616C 3 2752vr 4 3584 2 , 2 1024 ln 4 2 

H 7r m2 H 1 7r In 2 H 

81 9 27 135 9 9 

+ 8192a 4 \ J 4081770917 875224tt 2 199360^ ^ 1395520C 3 



3 J \ 279936 243 27 81 

97696tt 4 127232 9 , 2 36352 ln 4 2 290816a 4 11008 . 

+ — vr 2 ln 2 2 + + + — 7T 4 ln2 



405 27 27 9 45 

14336 2, 3a 4096 ln 5 2 nnrr/%a 1792tt 2 C 3 87296C 5 * 



9 



i 2 3 4096 ln 5 2 1792tt 2 C 3 8^96Cs\ 
-tt In 2 - + 32768a 5 - — —j 



- / 125873914573 13545868tt 2 3500896 2 24506272C 3 428624vr 4 

'' 1679616 + 729 81 ^ U + 243 243 

2791040 2l 2 797440 ln 4 2 6379520a 4 390784 . 508928 2 , , 

+ 7r ln 2 2 + + + vr 4 ln2 tt 2 ln 3 2 

81 81 27 135 27 

145408 ln 5 2 1163264a 5 63616tt 2 C 3 3099008Cs 93184tt 6 22016 7r 4 ln 2 2 



45 3 27 9 405 15 

14336 9l 4 8192 ln 6 2 7168 2l , 633344(1 

H — ix ln 4 2 + — + 393216a 6 + 180224s 6 + -— tt 2 ln2 Cs 7 --^ 

3 15 3 9 

6 / 3593750577461 193770934tt 2 54183472 2l _ o , 379284304C 3 



' ~ " 10077696 + 2187 243 ^ ln2 + 729 

37634632tt 4 49012544 9 , 2 14003584 ln 4 2 112028672a 4 1714496 . 

+ tt 2 ln 2 2 + + + tt 4 ln2 



3645 243 243 81 81 

11164160^ 2 ln 3 2 637952 ln 5 2 25518080a 5 1395520tt 2 C 3 67981760C 5 



81 ~ 27 9 81 27 

3308032tt 6 781568 4 , 2n 508928 2l 4rt 290816 ln 6 2 

-7T ln 2 2 + — - — 7r ln 4 2 + + 4653056a 6 



1215 45 ~ 9 45 

6397952s 6 254464 2l > 22483712(1 719872 fil 88064 4 , , 

+ + vr 2 n2 C 3 — + tt 6 n2 + vr 4 n 3 2 

3 9 S 27 405 15 

57344 2l , 32768 ln 7 2 19922944, n 1789696tt 4 C 3 
vr 2 ln 5 2 + 4718592a 7 ln2 s 6 + 



35 21 189 

1JW ,, 2n , 24903680, ^ >2 38040320tt 2 C 5 nnn „ n *,2 n + 72259840C 7 
-143 367T 2 ln 2 2 ( 3 + ln2 C 3 + — + 2095104 ln 2 2 ( 5 



21 63 7 

35C 

2] 21 



, 19922944s 7a 25493504^ 7 , 

+ — ™ ^ — +°( e n> ( 34 ) 
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- £ (-f + ^^ + ^)-,(- 734 -^ +20C 3 

551tt 4 76tt 2 C 3 a \ ,/ 12254 775tt 2 1232C 3 

+ + — + 462C 5 - e 3 

90 3 S 7 V 3 3 3 

28^ + 236^, + i482(s + 2353^ + 482^ _ ^_60346 

_ 1383 ^ - ^ _ ^ _ 32^ C 3 - 252 Cs + 

3 45 S S 1890 

+ i_5ioa + 357^3 + 8M7T% + 15307C ^ + (e5) I (35) 



.. , 2 10 i /26 vr 2 \ / IItt 2 16Cs~ 
G 5 , 2 (4 - 2e) = r(l + __ - _ _ e -^_ + _ j - (^2 + — + 

V 3 3 3 45 J V 

,™ 2, rt I888C3 3 vr 4 128 2l 2rt 641n 4 2 _ 8tt 2 C 3 
-160vr 2 ln2 + — - — + — vr 2 ln 2 2 + — - — + 512a 4 + — — + 96C 5 
3 5 3 3 3 

0/ 17470 1817tt 2 2l > 751tt 4 1280 2l 2 6401n 4 2 
-e 3 — + — 1024vr 2 Lq2 + 3600Cs tz- + ln 2 2 + 



45 3 3 

736 1024 2l o 5121n 5 2 8tt 2 C 3 , 368tt 6 

+5120a 4 + — vr 4 ln2 - — — 7T 2 ln 3 2 — + 4096a 5 - 2496C 5 + —r^ 

45 9 15 3 945 

64C 3 \ 4/ 85562 „ „ „, 53264C 3 5849tt 4 8192 2l 2 

+— ^ J - e 4 + 2649tt 2 - 5376tt 2 ln2 + + vr 2 ln 2 2 

3 J V 3 3 45 3 

4096 ln 4 2 „ 1472 . 10240 2l , 1024 ln 5 2 

+ + 32768a 4 + tt 4 ln2 tt 2 ln 3 2 + 40960a 5 

3 9 9 3 

3 76tt 2 C 3 2747T 6 3424 4 , 2144 , 4 2048 ln 6 2 

3— " 28512Cs " IT" " IT* ln2 +^ 1 2 n2+ ^45— 

+2567r 2 a 4 + 32768a 6 + 12288s 6 + 352tt 2 ln2 ( 3 - i-^^ 3 + l^lli 3 . + 4 87r 2 ( 5 

3 45 



+ 1328Cr) +0(e 5 ) I , 



(36) 
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G 7 (4 - 2e) = r(l + e) 3 | - 5C 5 ) + e(ln% - 10( 5 + ^ + 7C 3 2 

+e 2 - 20C 5 + ^ + 14 C 2 - + 166^5 - 212( 7 

+0 (e 3 ) } , (37) 



+ 200C3 + 37^ _ \ _ ^ /_ 1244 + igg7r2 + _ 46^ _ 2i7r2 ^ + 



n .n 2, « 80 2, 2n 6411142 r1n A 3^ „r„n 3100tT 2 ro „ ft/ . 218tT 4 

-168tt 2 ln2 + — 7r 2 ln 2 2 + — - — + 512a 4 - e 3 -7572 + — - — + 5360C 3 zr- 



10 ) \ L5 

80 2l 2 641n' 

— 7T 2 ln 2 2 + 

3 3 

>2tt ( 

+280tt 2 ln 2 2 - 6vr 4 ln 2 2 - ^vr 2 ln 3 2 + 224 ln 4 2 + 6vr 2 ln 4 2 - !H!^jL + 5376a4 



.^!<3 _ ig76C5 _ 22^ _ 332 ^ 2 _ ii2g ^ 2 ^ 2 + 128^ ln2 + ^^2 ln2 



+ 144vr 2 a 4 + 3072a 5 ^ + O (e 4 ) I , 



(38) 
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SM(4 _ ae) . IV + «,^ + ^ + ^(»_^) + (»_ a H_^_^ 

2 I6C3 4577T 4 267T 2 C 3 OAn , , 14717T 6 , 8 41 2o 

— 16n m2 H 342Cs H h -7r In 2 

3 90 3 S5 2835 3 

8 .2, 4 «/(„2„ RfiJ|„o/ , «OA2\ , a/ 26335 2 2 



-vr 2 ln 4 2 - 647r 2 a 4 - 56tt 2 ln2 Cs + 82tfj + e 3 ^— — + 187tt 2 - 224vr 2 ln2 

rno/ . 277vr 4 224 2l „ 641n 4 2 r _ 58tt 2 C 3 , fton/ . 

+598C 3 - — 3- + ^vr 2 ln 2 2 + — - — + 512a 4 ^ - 1082C 5 

Id o o o 

299tt 6 . 9 „, , 9 9 400G 2 112 fi 

+ 8tt 4 ln 2 2 - 8vr 2 ln 4 2 - 1927r 2 a 4 - 168vr 2 ln2 Cs + - -77 vr 6 ln2 



378 " " ~ -• - ... 135 

-^tt 4 ln 3 2 + ^tt 2 ln 5 2 + 768tt 2 ln2 a 4 + 768vr 2 a 5 + 1024 ln2 s 6 
3 5 

165l7r 4 C 3 1232 „, 9 > 224, . _ , >9 499tt 2 C 5 
+ 77 + ^7-^ 2 ln 2 2 Cs - ln 4 2 C 3 - 1792a 4 Cs - 1280 ln2 Cf ^ 

(39) 



-8365C 7 - 1024s 7a - 1024s 7 ^ + O (e 4 ) J , 



. 1 7 31 /103 7T 2 2C 3 4tt 4 \ /235 

G6 , ( 4-2 £ ) = r(i + £) ^- + - + - + (- + T + ^--j +f (- 

14tt 4 rt 2 , 10097T 6 8 4l o 8 2l . „, , „ n> 

+ + 2tt 2 Cs + 6C5 + j^ggg + 3 71 " ln 2 - 3^ ln 2 " 6471 ° 4 - 56?r ln2 & 

a2 \ ,/ 3953 2 ~ 1760Cs 307vr 4 224 2 , 2 

+32C? + e 3 + 186tt 2 - 224vr 2 ln2 + + + vr 2 ln 2 2 

S3 J V 3 3 90 3 

641n 4 2 9> , 979tt 6 , 2 9l 4 9 

+ + 512a 4 + 14vr 2 Cs + 222C 5 + + 8tt 4 ln 2 2 - 8tt 2 ln 4 2 - 1927r 2 a 4 

3 630 

-1687T 2 ln2 Cs + — - — tt 6 ln2 - -tt 4 ln 3 2 + — tt 2 ln 5 2 + 768vr 2 ln2 a 4 
3 135 3 5 

+7687r 2 a 5 + 1024 ln2 s 6 + — — + -— vr 2 ln 2 2 Cs - — ln 4 2 Cs - 1792a 4 Cs 

45 3 3 

-12801n2C 3 2 - 215tt 2 C 5 - 6114C 7 - 1024s 7a - 1024s 7f ^ + O (e 4 ) J , (40) 
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., , 1 3 t /55 tt 2 \ /95 7tt 2 14Cs 4tt 4 ^ 
G 6 , 3 (4-2 e )=r(l + e) 3 — + — + e- 1 --- + - 



6e 3 2e 2 V 6 3 / V 2 3 3 45 

/1351 31tt 2 „ 2l n 23vr 4 25tt 2 C 3 49( 5 \ 2 /2023 

^ + 32l2 , n2 _ _ 148^ _ 40 r2 ^ _ 201^2 _ . 109.% 



3 ~ 3 45 3 3 6 

413Cs 237r 6 7 4l „ 7 2l , r „ , , ft 2l n/ . 191CI 
+ — — + -vr 4 ln 2 2 - -vr 2 ln 4 2 - 56vT 2 a 4 - 49vr 2 ln2 Cs + 



2 54 3 '"' ~ 3 ' ~ 4 

,/26335 235tt 2 _ 371tt 4 2l 2 , 4 

+e 3 + 140/t 2 Lq2 - 994C 3 80tt 2 ln 2 2 - 56 ln 4 2 

\ 6 3 30 

,o,, 182 4i « 368 2, 3a 1841n 5 2 1297r 2 C 3 107C 5 

-1344a 4 - — vr 4 Lq2 + — vr 2 ln 3 2 + 1472a 5 + — ^ + — ^ 

4o y lo ii 

457tt 6 35 4l 2rt 35 2l 4rt rtnrt 2 rtJr 2l _ . 2641(1 
+ W + T ~ ¥ ~ 4 ~ C3 + ~12 

-^tt 6 ln2 - ^tt 4 ln 3 2 + ^tt 2 ln 5 2 + 656vr 2 ln2 a 4 + 656vT 2 a 5 + 800 ln2 s 6 
54 9 15 

3803vt 4 C 3 , 1036 _ 2 2oA 175, 4oA A innnl _ n/ .2 2 057vt 2 C 5 



+ 1Qn + —7^ ln 2 2 C 3 - — ln 4 2 Cs - 1400a 4 Cs - 1000 ln2 C| - 
800s 7a - 800s 76 



180 3 3 ^ ^ 12 

44413C 7 



8 



>)+0( e 4 )j, (41) 



G M (4 - 


2e) = r 


41vT 4 


^Cs 


90 


3 


+42Cs - 


1447tt 6 


2835 



„ \ 2 / 119tt 2 _ 2l rt ooA 14vr 4 26tt 2 C 3 
+ 44C 5 ) - e 2 f — + 16vr 2 ln2 + 38C 3 + — + — ^ 

- -vr 4 ln 2 2 + -vr 2 ln 4 2 + 64vT 2 a 4 + 56tt 2 ln2 Cs - 54C 3 2 
3 3 

o/796vr 2 nni 2l n 51tt 4 224 2l „ 641n 4 2 „ n 

+e 3 224tt 2 ln2 + 256C 3 + + tt 2 ln 2 2 + + 512a 4 

V 3 10 3 3 

31vT 6 

-26tt 2 C 3 - 574( 5 + + 8tt 4 ln 2 2 - 8tt 2 ln 4 2 - 192vT 2 a 4 - 168tt 2 ln2 Cs 

630 

112 64 1 28 

+104Cf 7T 6 ln2 vr 4 ln 3 2 + tt 2 ln 5 2 + 768vt 2 ln2 a 4 + 768vT 2 a 5 

1 35 3 5 

93l7r 4 C 3 1232 „, 2 > 224, 4 > 
+1024 ln2 s 6 + — — ^ + -^n 2 ln 2 2 Cs - — ln 4 2 Cs - 1792a 4 Cs 
45 



-1280 ln2 C 3 2 - 217vt 2 C 5 - 5655C 7 - 1024s 7a - 1024s 76 ^ + O (e 4 ) J , 



(42) 
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., , 1 7 31 /103 4( 3 2vr 4 \ /235 
G 6 , 5 (4-2 e )=r(l + e) 3 — + _ + -+ _ +ef 



3e 3 3e 2 3e V 3 3 15 J \ 3 

8tt 2 32C 3 3tt 4 28tt 2 C 3 „ oa \ 2 /19 112tt 2 nrt 2l rt 692C 3 
+ + + — - 78C 5 + e 2 — + 32vr 2 ln2 + 



335 3 '7 V 33 3 

164tt 4 16 2l 2 161n 4 2 140tt 2 C 3 , ^ 928tt 6 
- — 7T 2 ln 2 2 + — - — + 128a 4 + r-^ - 414C 5 + 



45 3 3 3 ^ 945 

+— vr 4 ln 2 2 - — 7T 2 ln 4 2 - 1287T 2 a 4 - 112tt 2 ln2 C 3 + 169C| + e 3 

o o / 

-4487T 2 ln2 + - ^ + 967T 2 ln 2 2 + 96 ln 4 2 + 2304o 4 + ^tt 4 ln2 

32 2l o 321n 5 2 w „ n 532vr 2 C 3 946vr 6 80 4l 2rt 

+ tt 2 l n 3 2 + 768a 5 + ^ - 2246Cs + —— + — ^ ln 2 2 

o O o lay o 

80 o, 4^ 2 ^« 2, 2519CI 266 fil rt 128 , 

-— vr 2 ln 4 2 - 6407r 2 a 4 - 560vr 2 ln2 ( 3 + — ^ - — tt 6 La2 - — tt 4 ln 3 2 
3 3 135 3 

i n s 2 + 15367T 2 ln2 a 4 + 15367r 2 a 5 + 2432 ln2 s 6 + + — vr 2 ln 2 2 & 

5 5 3 



ln 4 2 C 3 - 4256a 4 C 3 - 3040 ln2 C 3 2 - 132tt 2 C 5 - _ 2432s 7a - 2432s 7fe j 

+0 (e 4 ) } , (43) 
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G 6 , 6 (4 - 2e) = r(l + efi ^ - {% - 10& + ^ + 1„2 



90 } V 3 

13tt 4 49tt 2 C 3 85Cs\ 2/^2 „ 21 n ,0, 2177r4 
+24C 3 + — — - e 2 -7vr 2 - 16tt 2 ln2 - 42( 3 + 



18 6 2 / V 90 

,40 2 2o 201n 4 2 245tt 2 C 3 , 425( 5 1751tt 6 4 2 
+y vr ln2+^ + 160a 4 - — — + ~ 2 18f ~ 5n ln 2 

r 2 ln 4 2 + 1207r 2 a 4 + 105tt 2 ln2 6 - + e 3 ( - 84vr 2 



+5^ ln 4 2 + 1207r^a 4 + 105^ ln2 C 3 + e d ^— 84tt" ln2 

223vr 4 16 2l 2o 1041n 4 2 ooo 182 . n 368 2l 3rt 

+268C 3 vr 2 ln 2 2 832a 4 tt 4 ln2 + tt 2 ln 3 2 

S3 18 3 3 45 9 

, 1841n 5 2 313tt 2 C 3 431( 5 , 1751tt 6 4 2 „, 2| 4 „ 

H 77 1472a 5 H 1 1 — h 25tt ln 2 - 257r ln 2 

lo z z 610 

-6007TV - 5257T 2 In2 (3 + - ^ 6 ln2 - ^tt 4 ln 3 2 + ^vr 2 ln 5 2 

4 90 9 15 

+1424vr 2 ln2 a 4 + 14247r 2 a 5 + 1824 ln2 s 6 + 925971 ^ 3 + 756tt 2 ln 2 2 ( 3 - 133 ln 4 2 ( 3 

180 

, >9 3841tt 2 C 5 105133C 7 

-3192a 4 Cs - 2280 ln2 C 3 2 — — — — - 1824s 7a - 1824s 7fe 

12 8 

+0 (e 4 ) 1 , (44) 
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G 8 (4 - 2c) = r(l + ef \ [ ~ + 4vr 2 ln 2 2 ) + e ^ - 16tt 2 ln 2 2 + ^vr 4 ln2 

208 9 , , 321n 5 2 I7n 2 ( 3 \ 9 f 14vr 4 „, 2 

vr 2 ln 3 2 + 256a 5 + — - 291C 5 + e 2 + 112tt 2 ln 2 2 

9 15 5 6 S 7 V 3 

664 . o 832 2l 3o 1281n 5 2 innA 34vr 2 C 3 

-vr 4 ln2 H 7T 2 ln 3 2 H 1024a 5 — + 1164C 5 



45 9 15 3 

21743vr 6 197 2 713 , 641n 6 2 „ 

— tt-tz —7i 4 ln 2 2 + --7r 2 ln 4 2 + — 104vr 2 a 4 + 5120a 6 + 2688s 6 

1134U 9 9 9 

9 , >A ,/80tt 4 a „, , 4648 4l 5824 _ , 

-51vr 2 ln2 Cs - 953C, 2 + e 3 640tt 2 ln 2 2 + vr 4 ln2 vr 2 ln 3 2 

3 7 V 3 45 9 

8961n 5 2 m 238vr 2 C 3 A 21743tt 6 788 4l 2 2852 2n 4 

+ 7168a 5 + — - 8148C 5 + + vr 4 n 2 2 vr 2 n 4 2 

15 3 S 2835 9 9 

~ >b hl 2 + 416vr 2 a 4 - 20480a 6 - 10752s 6 + 204tt 2 ln2 Cs + 3812C 3 2 + -^vr 6 ln2 



9 ' 189 

+ ln 3 2 _ 7288 2 5 _ 4864 ln 7 2 2 i 5 2 07 r 2 a 5 + 77824a 7 

135 45 315 457 

106880, 40037r 4 C 3 2167 2l 2 , 316 , 4 _ 
— ln2 s 6 + ^-^ + —tt 2 ln 2 2 ( 3 - — ln 4 2 ( 3 - 2528a 4 C 3 

133600 9 . 2 ,8755617^5 2 9 , 1325727C 7 , 106880s 7a 

H — ln2 C 3 H — h 37200 In 2 ( 5 1 

7 3 84 7 7 

161920s 76 \ 



(45) 



Here 



a n = Li n (l/2) , 

se = C-s.-i + Ce , 

s 7a = C-5,1,1 + C— 6,1 + C-5,2 + C-7 

•57b = C7 + (5,2 + C-6,-1 + C5-1.-1 



4 Conclusion 

Using the DRA method [1] we have performed the analytic evaluation of two families 
of the three-loop master integrals in the e-expansion up to the transcendentality 
weight intrinsic to the corresponding four-loop master integrals. This method is 
naturally combined with other methods. First, it hardly relies on IBP reduction 
which is necessary in order to obtain difference equations with respect to dimension 
in a closed form. To do this we used the code based on Ref. [24J and the code called 
FIRE [22]. To reveal the position and the order of the poles in a basic stripe we used 
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a sector decomposition [261 123 [28] implemented in the code FIESTA [281 122]. To fix 
remaining constants in the homogenous solution of dimensional recurrence relations 
we applied the method of Mellin-Barnes representation [3Ql EH [32] . When dealing 
with multiple zeta values we used the code HPL [33]. Finally, to arrive at analytical 
results for coefficients in the e-expansion we applied the PSLQ algorithm [20] . 
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